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Abstract 

We calculate the kaonic contribution to the shear viscosity of quark matter in the kaon-condensed 
color- flavor locked phase (CFL-KO). This contribution comes from a light pseudo-Goldstone boson 
which arises from the spontaneous breaking of the flavor symmetry by the kaon condensate. The 
other contribution, from the exactly massless superfluid "phonon" , has been calculated previously. 
We specialize to a particular form of the interaction lagrangian, parameterized by a single coupling. 
We find that if we make reasonable guesses for the values of the parameters of the effective theory, 
the kaons have a much smaller shear viscosity than the superfluid phonons, but also a much shorter 
mean free path, so they could easily provide the dominant contribution to the shear viscosity of 
CFL-KO quark matter in a neutron star in the temperature range 0.01 to 1 MeV (10*^ to 10^*^ K). 



I. INTRODUCTION 



In this paper, we explore the shear viscosity of one of the predicted phases of high 
density quark matter. Transport properties of quark matter, such as the shear viscosity, 
are of interest because they are the basis for signatures by which we could infer or rule 
out the presence of exotic phases in the core of neutron stars. Previous work on transport 
properties has studied the bulk viscosity [1-4], specific heat [5, 6], neutrino emissivity [7-11], 
and thermal conductivity [12-16]. 

The shear viscosity is phenomenologically relevant because it damps physically important 
excitation modes of the star. In particular, a fast-rotating neutron star will spin down rapidly 
if the internal viscosity is too low, because of the spontaneous growth of r-modes [17]. The 
observation of rapidly rotating neutron stars can therefore be used to place limits on the 
internal viscosity. Calculating the viscosity of candidate phases then allows us say whether 
these phases can be present in the neutron star. 

This paper studies the kaon-condensed color-flavor locked phase (CFL-KO) of quark mat- 
ter. This is a candidate phase of quark matter at the central density of a neutron star. 
For a review, see Ref. [18]. The CFL-KO phase has a condensate of kaons [19, 20], which 
spontaneously breaks the flavor symmetry, producing a Goldstone boson. Since the flavor 
symmetry is also explicitly broken by the weak interaction, this Goldstone kaon acquires a 
small mass in the keV range [21]. Our analysis is relevant to temperatures above this value, 
where one can ignore this small mass. We use the effective theory of the Goldstone kaon, 
which was worked out in Ref. [5]. The full interaction lagrangian has three independent 
coupling constants, but we will specialize to a specific ratio of their values, which makes 
our results dependent on one overall kaon interaction coupling. This enables us to make an 
estimate of the expected scale of the shear viscosity in this phase. We defer the calculation 
for the most general interaction lagrangian to future work. 

As well as the shear viscosity we calculate the mean free path of the Goldstone kaons, 
since the concept of shear viscosity will only be applicable to the matter in the neutron star 
if there is local equilibration on distance scales much smaller than the size of the star. 

The rest of the paper is laid out as follows: Section II will discuss the low energy effective 
theory of the CFL-KO phase and the interactions among the lowest energy excitations of 
the theory. Section III will show the results for the contribution of the kaons to the shear 



1 



viscosity. Section V will present some conclusions and discuss future directions of this work. 
In the appendices we cover technical details of the calculation of the mean free path, the 
treatment of the co-linear regime of the collision integral, and the approximate evaluation 
of the collision integral. 



II. LOW ENERGY EFFECTIVE THEORY 



A. Lowest-order lagrangian 

The low energy degrees of freedom in color-fiavor-locked phases of quark matter are the 
massless superfluid Goldstone mode, arising from the spontaneous breaking of baryon num- 
ber, and the light pseudo-Goldstone meson octet, arising from the spontaneous breaking of 
three-flavor chiral symmetry. The contribution of the superfluid mode to transport proper- 
ties has been studied previously [3, 22]. We focus on the contribution from the meson octet, 
described by a meson field S whose effective lagrangian up to second order is [19, 20] 

L= ^TrfDoSDoS^ -w^VSVE^] +a4^detMTr[M-^(S + S^)] (1) 

where DqL = doT, — i[A,T]. The Bedaque-Shafer effective chemical potential [19] is A = 
— ^j^, l^q is the quark chemical potential, and M = diag(m,j, m^, m^) is the quark mass 
matrix. At asymptotically high density the constants /tt, v, and a can be determined by 
matching the effective theory to perturbative QCD, [23, 24] 

21 - 81n2 /i^ , 1 3A2 

where A is the fermionic energy gap at zero temperature. This dependence of a on A is also 
seen in NJL models [25], so from now on we will work in terms of /,r and A, assuming that 
a is given by (2). The meson field T, can be parameterized in terms of fields 6a, 

S = exp{z9/U) , (3) 

where 6 = OgTa-, and are the Gell-Mann matrices of SU(3) with normalization tr(TaTb) = 
25ab- The and are the lightest mesonic degrees of freedom [23, 24], and electric 
neutrality disfavors the presence of charged kaons (since they must be balanced by electrons). 
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so we focus on the neutral kaons, and K^, corresponding to 6q and 6j. The zero- 
temperature neutral kaon mass and chemical potential can be deduced from the Lagrangian 



m]^ = amuirrici + m 
ml-m^ 



2 _2 (4) 



We will assume that > i^k, so there is kaon condensation. The critical temperature for 
kaon condensation is expected to be of order tens of MeV, well above typical neutron star 
temperatures, and of the same order as the critical temperature for the CFL condensate itself 
[5]. We will also assume, following Ref. [5], that the condensate is small, so fiK is only a little 
larger than rriK- Given the uncertainty in the effective theory couplings, this is a perfectly 
conceivable scenario: taking for example fi = 400 MeV, = 100 MeV, A = 80 MeV and 
rriu = 5 MeV (and ignoring because it makes a negligible contribution) one obtains 
fiK = 12.5 MeV and rriK = 11.0 MeV. It is then convenient to define, following Ref. [5], an 
energy gap 

I Stti I 

6m = rriK — , ^ 1 . (5) 

Note that 6m is negative in the CFL-KO phase. Because \6m\ ^ 1 we can usually treat 
and mx as being identical to leading order in 6m (an exception is discussed in Sec. IV). 

A self-consistent calculation [5] (see also [26]) then yields the excitation energies in the 
neutral kaon sector, 



Ei = + /i^ T yJ^f^lE^p + SM^ , (6) 

where 

= vY + W . (7) 

In that self-consistent calculation, M and 6M were thermal masses that depended on tem- 
perature and the underlying mass and chemical potential (see Eq. (81) in Ref. [5]). In this 
paper we are interested in low temperature range applicable for compact stars. In this case, 
the thermal masses become independent of temperature and are given by 

= 2/i|- -m]^^m\ , 
6M'^ = iJ.K-m%^ 2mK\6m\ . (8) 

The mode with energy E_^ is massless: this is the Goldstone kaon. We can define a corre- 
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spending field using the parameterization 

9q{x) = (0 + p(x)) COS'i9(x) , 

e^{x) = {(f) + p{x)) sini9(x) , (9) 
i^ix) = /^sin(0//^)^9(x) , 

so (j) is the kaon condensate, p is the massive radial mode, and is the angular Goldstone 
mode which we have then rescaled to make a scalar field ip with a canonically normalized 
quadratic derivative term and the conventional energy dimension of 1. The mass of the 
radial modes is given by the value of E^{p = 0) = ^/QjJ^j^ — 2m^- ^ 2mK, which is typically 
on the order of a few MeV, so at the 10 to 100 keV energy scale, which is relevant to neutron 
stars, it is heavily suppressed and can be ignored. The magnitude of the kaon condensate is 
[5] 



^).4/.^M (10) 

We can then linearize Eq. (6) to obtain a linear dispersion relation for the Goldstone 
kaon, 

E{p) = up 



The error involved in approximating Eq. (6) by Eq. (11) is less than 5% for 
p < 0.6 mK\Sm\/v. 

We will obtain the contribution to the shear viscosity from the Goldstone kaon. For 
this we need its interaction lagrangian, but it is easy to see that (1) does not contain any 
interaction terms for the field ip. This follows from the fact that ip, as a Goldstone boson, 
must couple via derivatives, and (1) only goes to second order in derivatives. We therefore 
need to write down higher order derivative terms in the effective theory to obtain interactions 
among the Goldstone modes. 



B. Interaction lagrangian for the Goldstone kaons 

We obtain higher derivative terms in ip by writing down the leading higher derivative 
terms in the lagrangian for S, and using (9). We keep only terms with the symmetries of 
the system, namely rotational symmetry, parity, time-reversal, and the SU{3)l ® SU{3)r 
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la. (TV[DoSDoSt])' 



lb. 2TT[{DoJ:DoJ:^f 



J TT J TV 



2a. (Tr[VSVEt])^ ^ 
2b. 2Tr[(Vi;VSt)2] 

3a. TY[L)oSDoSt]TV[VSVSt] 4,^ ,,2,^, ,2 8/i/<sm(^ 



3b. 2Tr[L>oSDoStVSVi;^] 



TABLE I: The six interaction terms at fourth order in derivatives for the effective theory (first 
column), and the interaction terms for ip that they transform to using (9), when terms involving p 
are dropped. In the effective lagrangian they have coefficients of order /^/A^. 

chiral flavor symmetry. We also discard terms that, when we substitute (9), will produce 
interactions that all involve the p field; an example is three- derivative terms where S enters 
four times. The allowed terms with no more than four derivatives of S are shown in Table I 
(left column). Since the effective theory breaks down at momenta of order A (for example, 
scattering of Goldstone bosons at that momentum will produce quasiquarks, which are not 
included in the effective theory) we expect that the momentum expansion will be in powers 
of (1/A)V [19]. We therefore expect the interactions in the left column of Table 1 to occur in 
the lagrangian with coefficients Cj/^/A^, where the Cj are dimensionless coupling constants. 

Using (9) and dropping terms that involve the heavy field p, these six terms reduce to 
the three corresponding interaction terms for ip shown in the right column. In each case, we 
find two different interaction terms for S reduce to the same interaction term for ip. This 
means that the interaction lagrangian for ip only depends on three linear combinations of 
couplings. Note that in Table I we have defined a scaled version of the kaon condensate 
expectation value (9), 




^^<P/U^2\^ (12) 
rriK 

The interaction Lagrangian for ip can then be written out as 

^ = C^i;^(5oV^)^ + C^3-^(9oV^)^W)^ + C72-^(VV^)^ 



(13) 



At this point we specialize to a particular form of the interaction lagrangian, with the 



following relationship among the three coupling constants. 

C = C, = C2 = -\C,. (14) 

This reduces the number of coupling constants from three to one. The remainder of our cal- 
culation is for this special case, adopted because it leads to a particularly simple interaction 
lagrangian which is similar to that written down for the superfluid phonon in Refs. [22, 27], 



Ant = -^,{d,i,d^^ljf + -^(9oV^)(a^^a'^^/^) (15) 



where 

A = 8C^ and ^ = 16 sin((^)C (16) 



We leave the analysis of the fully general interaction lagrangian (13) for future work. 



III. SHEAR VISCOSITY 



The shear viscosity 77 is a coefficient in the viscous stress tensor 5Tij , which is the deviation 
from equilibrium of the momentum flux tensor Tjj for a fluid with pressure P and energy 
density e, 

= -r]V,, + ■■■ (17) 

where 

V,j = Vj + d, V, - '^6ijV ■ V (18) 

and the ellipsis in the equation for STij stands for other dissipative terms arising from 
phenomena such as bulk viscosity and thermal conductivity. V(x, t) is the fluid velocity at 
a given position and time. We will need to make sure that we only need to keep the first 
order dissipative terms in Eq. (17) (see discussion at the end of IV B). The stress-energy 
tensor and the viscosity can be calculated using kinetic theory [28] . For a system of identical 
bosonic particles with dispersion relation Ep, 
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where v is the velocity of the Goldstone kaon (see Eq. 11). The full distribution function is 
given by 

^) = ,....(x.!)/T _ 1 = n + t) (20) 

where m'^(x, t) is the 4- velocity of the fluid, and Sfp is a small departure from the equilibrium 
Bose-Einstein distribution 

For shear viscosity we are interested in deviations from equilibrium arising from a shear 
flow, so we write the deviation from equilibrium as 

SUM = -^^^f,{i + 

X{p,^,t) = g{p)pMVki{x,t), (22) 

where 

Pki=PkPi-^Skip'^. (23) 
Substituting (22) into (19) and (17) we find 



5r.,(x,t) = -u^ J ^/-^'^^m + /°) M^,t). (24) 

Using the definition of Vij (see Eq. 18) one can write 6Tij (Eq. 17) in the following form 

r] 2 
STij = --[SikSji + 6ii6jk - -Sij6ki]Vki, (25) 

Comparing this to Eq. (24) gives us 

^TA A A 2. . . of _fPPiPj9(j>)PklrQ(. , .0^ /^^x 

-[dikdji + dudjk - g^ijM = ^ y (2^~E f — ^p'' ^ ^ 

Then, by contracting the tensor on the left hand side with respect to the pairs of indices i, k 
and j, / we can determine the shear viscosity in terms of the function g{p), 

^=1^1 P'fvi^ + OaiPl (27) 



15 T 

where we have adopted the notation 



2Ep{2nf' ^^^^ 
Using the fact that = ^PijPij (see Eq. 23) one can write an alternate form of the shear 
viscosity which will be used later, 

21.2 



5T 



/ fpil + fp)9ip)P^,P^J■ (29) 
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To solve for the viscosity, we need to find a form for g{p). To do so, we can use the 
Boltzmann equation given in the absence of external forces by 

f ,V.V/,^C|/J. (30) 



The left-hand side can be written as [28] 

dfp f p 



[^ + OPijVij . (31) 



dt 2pT 

Again, this specific form assumes that we are only interested in shear flows. Thermal 
gradients and bulk flows would give additional term on the right hand side of Eq. (31). It 
is this form that also helped motivate the structure of the ansatz in Eq. (22). 

The collision operator C[fp] should contain any possible collision terms for the kaons. We 
will restrict ourselves to the terms lowest order in the coupling constants as more vertices are 
suppressed because each vertex brings in more powers of 1/ fn or 1/A (see (15) and (16)). 
Also, we will ignore the 1^2 processes because for a particle with a linear dispersion 
relation such processes must be co-linear, so they do not involve momentum transfer that 
would contribute to the shear viscosity. Finally we are left with the collision operator for 
2-body scattering given by [27] 

C2^2 = 7^ [ {2Tif5\P + K-P'- K')\M\^D^^^ (32) 

"^^p Jk,k',p' 

Where Ai is the 2 — )■ 2 scattering amplitude and -02^2 contains the distribution functions. 

D2^2 = fp' fk'{l + /p)(l + fk) - fpfk{l + fp'){l + fk'). (33) 

We can also linearize the distributions as D ~ D^ + 6D using our definition of Sfp in Eq. (22). 
-D° would make the collision integral vanish by detailed balance. One can then write 

SD2^2 = ^f^mi + /;0(1 + fk') ixip) + X{k) - Xip) - X{k')) (34) 
and the collision integral as 

C2^2 ^ S^l {27rmP + K-P'- K')\M\'f,{l + /°)(1 + f) 

^^p-L Jk,p',k' 
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Using (31) and the Boltzmann equation, we can conclude that 

'^'fA^ + Op^^ = F^A9{p)] (36) 

One can then use this equation and (29) to get another expression for the shear viscosity in 
terms of coUision term 

^ = ^ ^ 9{p)PijFij [gip)] (37) 

The process now is to evaluate rj from (27) and (37), and ensure that they give the same 
answer. Formally, this is equivalent to solving the Boltzmann equation directly. Ensuring 
that the two forms are equal is quite non-trivial and is typically done by expanding g{p) 
using an orthogonal set of functions, [29-31] 

N 

9{p)=p''Y.^sBs{p) . (38) 

This expansion introduces two new parameters: A^, the order of the polynomial approxima- 
tion; and n, which we call the minimum-exponent parameter, because the lowest power of 
p that occurs in the polynomial expansion is p". The correct result is obtained in the limit 
N oo for any value of n. However, as we will see, the rate of convergence with increasing 
is strongly dependent on the minimum-exponent parameter n. 

The polynomials Bs{p) are defined such that the coefficient of the highest power is 1 
and they obey the orthogonality condition [27] 



jp 



+ fl) Pr, P^JP''Br{p)B,{p) = ^^5,, . (39) 



These conditions uniquely specify the Bs{p) for all s, starting with Bo{p) = 1. From the 
orthogonality condition we find 

4"^ = 3 I + /;)/^" = ^;^r(6 + n)C(5 + n) . (40) 
Using g{p) from Eq. (38) in Eq. (27), and using the definition of aI^^ from Eq. (39), we get 

V = %boAt' . (41) 
An alternative expression for r/ follows from substituting g{j>) from Eq. (38) into Eq. (37), 

Af=oo 
s,t=0 



where 

Mst = ^ [ (2vr) V(P + K-P'- K')\M\'f;f,{l + /° )(1 + f}) p"P.(pKA*^. , 

^-^ Jp,k,k',p' 

= T^ j (2vr)V(P + K-P'- K')\M\'flfl{l + /° )(1 + §) A^^.A*^. 

-'-U-' Jp,k,k',p' 

(43) 

and 

A*^. = P*(p)p>., + Bt{k)k''hj - B,{k')k'%^ - Pt(p')P% • (44) 

The second hne of Eq. (43) uses the symmetry under relabehng the momenta of the legs in 
the scattering diagrams (P — K, etc.), and can be used to demonstrate that the diagonal 
elements of Mgt are positive definite. As we will see below, this ensures that the shear 
viscosity is also positive. 

Requiring the two forms of rj to be equal leads to a matrix equation for all the fej's. From 
that we extract b^, 

6o = |^4"HM"1)oo (45) 

where (M~^)oo means the first entry in the matrix inverse of Mgt. Using this in Eq. (41), 
we find 

V=^M^^?{M-')oo. (46) 

As noted above, this expression becomes accurate in the limit — t- oo, where the matrix 
M is of infinite size. It is known [22, 27, 29] that the result for finite A^ rises with A^, so for 
a matrix M^, with finite dimension A^, that obeys (43), 

V > |^(4"^)^(M^^)oo (47) 

We will see below that this expression converges rapidly with A^ for the optimal choice of 
the minimum-exponent parameter n. 

The remaining task is to evaluate the integral in Eq. (43). This requires the matrix 
elements for the 2^2 scattering process, iAi = iM.c + iM.s + i-Mt + i-Mu, see Fig. 1, with 
the individual channels being given by 

zMc = yA(P- K){P' ■K) + {P- K'){P' ■K) + {P- P'){K ■ K')] 

J n 

iMs = [2(po + ko)P ■ K + p,K^ + fcoP'] [2{p', + k'^)P' ■ K' + p'^K'^ + k'^P'^] G{P + K) 
iMt = iMs{P^K') 

iMu = iMt{P ^ K), (48) 

10 



FIG. 1: Feynman diagrams for the 2-body scattering amplitude. The black square is the 4-point 
contact interaction, the black circles are the 3-point vertices. 

where 

^(^) = T-2 2 2^ ^T m T (^9) 

is the Goldstone kaon propagator and the last two lines in (48) come from crossing symme- 
tries. 

The 12-dimensional integral in Eq. (43) can be simplified by eliminating the p' integral 
using the momentum-conserving delta-function. Then one can use the energy-conserving 
delta-function to eliminate the integral over the magnitude of k' . Three of the remaining 
8 integrals can be eliminated by selecting the z-axis to lie along the vector p, and noting 
that only the difference in the two remaining azimuthal angles matters. This leaves a 5- 
dimensional integral over the magnitudes of p and two polar angles corresponding to k 
and k' and one azimuthal angle. This can be evaluated numerically (see appendix C) using 
the Vegas Monte Carlo algorithm [32, 33]. 

The results that we present below are obtained by setting the minimum-exponent pa- 
rameter n to —1. This value is expected to give optimal convergence of the calculated shear 
viscosity to its physical value as a function of because, as shown in appendix B, this term 
most strongly suppresses the co-linear scattering and therefore give the smallest collision 
term. The shear viscosity is inversely proportional to the collision term and since we have a 
variational procedure that says the answer we get is a lower bound, we are only interested 
in the largest value of the shear viscosity that we can calculate. 

To check this reasoning we show in Table II results of calculations of the shear viscosity 
for different values of n and A^. We see that for n = —1 the value of rj at low is already 
close to the maximum (asymptotic) value at = oo. For n = —2 the convergence is almost 
as good, achieving ~ 1% accuracy at A^ = 2. For other values of n the convergence is 
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n 


iV = 


= 1 


iV = 2 


-3 


0.059 


0.861 


2.22 


-2 


1.59 X 10^ 


2.01 X 10^ 


2.02 X 10^ 


-1 


2.04 X 10^ 


2.04 X 10^ 


2.04 X 10^ 





0.610 


4.70 


19.4 



TABLE II: Values of shear viscosity in (MeV)"^ as a function of the order N of the polynomial 
approximation to g{p), for different choices of the minimum-exponent parameter n. The calculated 
value rises towards the physical result as — ?• oo, and in this limit should be independent of n. 
We see that for n = — 1 the result converges very rapidly as N rises, but for other values of n the 
convergence is slower. 

dramatically poorer. This behavior was also seen in Refs. [22, 27]. We conclude that we 
can achieve accuracy of better than 1% by choosing n = —1, and only keeping the first 
polynomial (A^ = 0), i.e. we set g{p) = 1/p. 

IV. RESULTS 

A. Analytic results 

We now describe how the shear viscosity depends on the temperature and on the param- 
eters of the effective lagrangian for the Goldstone kaon. The relevant parameters are the 
speed u (11) of the kaon and its interaction couplings A and g (15), which in turn depend on 
more basic parameters C, /^r. A, mx (16). Recall that we have assumed |5m| <^ 1, so in the 
expressions below, fix and rriK are usually interchangeable. One exception is the phonon 
speed which occurs in the shear viscosity raised to the lltli power (see discussion after 
(54)) so we use the full expression (the identity in (11)) for it. 

Before performing any numerical calculations, we can extract the temperature dependence 
of the shear viscosity. Because the co-linear scattering will not contribute to the answer, 
the propagator does not need to be regulated by the self-energy. Since the temperature 
only appears in the distribution functions and the self-energy, we can now factorize out the 
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temperature dependence by rescaling all the momenta by the temperature. Doing so we find 

where s and t are the indices of the matrix indicating how many terms we are keeping in 
our expansion for g{p) and n is the minimum-exponent parameter (see (38)). We also recall 
the temperature dependence of A^"'* from Eq. (39), 

Therefore, from Eq. (46), we obtain the temperature dependence of the shear viscosity 

cx T-^ (52) 

This dependence on temperature is also seen in other systems where viscosity arises from 
2^2 scattering of Goldstone bosons [22, 27]. The constant of proportionality in (52) has 
mass dimension 8. In the case of the shear viscosity due to phonons there was only one 
possible scale, the quark chemical potential /ig, so rju oc yU^/T^. However, we have several 
scales (/^, A, mx) manifesting themselves in two coupling constants A and g (16). Which of 
these is most important depends on whether the scattering is dominated by the contact term 
or by the exchange of a virtual particle. The dimensionless parameter u that determines 
which scattering process dominates is 

« = ^ = 96Csin^(^)(^)', (53) 



where (p = 1^\bm\lmK (12) and the 3 represents the 3 channels for virtual particle exchange. 
For typical values of (5m, m/^, and A, this ratio can be bigger or smaller than 1. When u ^ 1, 
the contact term dominates, so the scattering amplitude is proportional to A. When m ^ 1, 
the particle-exchange process dominates, so the scattering amplitude is proportional to . 
The shear viscosity is inversely proportional to the scattering cross-section, so 

1 /^^A^ 



M ^ 1 : r] = hi{u) 



M ^ 1 : r] = h2{iy] 



C2 



(54) 



where hi and h2 are dimensionless functions that depend only on the Goldstone kaon speed 
u, which depends on 6m and rriK- We can obtain their analytic form when u <^ 1. In that 
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case, the leading-order behavior in both regimes is a z/^^ power law (see appendix C). In 
general they must be calculated numerically, and the result (with the u^^ power law scaled 
out) is shown in Fig. 4. 

Finally, we note that the shear viscosity due to Goldstone kaons will be smaller than that 
due to the superfluid phonons, since /tt. A, and niK are much less than fig. 

B. Numerical results 

To begin, we will confirm the temperature dependence predicted in the previous section. 
To do so, we will fix the mass of the kaon, m^^ = 4 MeV, = 100 MeV, A = 100 MeV and 
C = 1. In Fig. 2, we show the viscosity as a function of temperature for a few values of 6m. 
The data points are obtained by numerical evaluation of the 5-dimensional integral, whereas 
the lines show the fit to a T^^ power law (52), which is independent of the regime of coupling 
constants' values. On the same plot, we show the contribution from the phonons [22]. As 
expected, the phonon shear viscosity is much larger. Most of the difference comes from the 
difference in magnitude of the coupling constants (the kaons coupling constant is larger) 
and the rest comes from a difference in the speed of the kaon and the phonon (the kaon's is 
smaller), which enters the expression for shear viscosity raised to a high power. Using the 
shear mean free path criterion p^^^^ < 1 km (appendix A) for the validity of hydrodynamics 
for neutron star oscillations, we expect the phonons to be non-hydrodynamic at T < 1 MeV; 
with the parameter values given above, the Goldstone kaons become non-hydrodynamic at 
T < 0.03 MeV. 

In Fig. 3 we show the shear viscosity as a function of the gap A, with Sm = —0.5 MeV 
and T = 1 MeV. The other parameters have the same values as in Fig. 2. This illustrates 
the transition between the two regimes given in Eq. (54). The crossover occurs at m = 1 
which corresponds to A = 30 MeV, which is indicated on the graph. As expected, we see 
that for large A (m <^ 1), 77 oc A^; for small A (m ^ 1), 77 oc A^. 

In Fig. 4 we present the results of numerical calculation of hi 2(1') (54). We have divided 
out the dominant behavior u^^ power law behavior (for details see appendix C). We see that 
the remaining z/ dependence is very mild, so to a good approximation the shear viscosity is 
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,24 








0.1 

T (MeV) 



1 



FIG. 2: (Color online) The shear viscosity as a function of temperature for kaons and phonons. For 
parameter values, see text. In the lower part of the graph, the points are numerical calculations 
and the straight lines are fits to the power law form given in (52). The phonons' calculated shear 
viscosity is many orders of magnitude larger, although using the shear mean free path criterion 
(appendix A) we expect them to be non-hydrodynamic in neutron stars at T < 1 MeV. 

given by (54) with 



To show how large the shear viscosity of CFL-KO quark matter could be, we look at a 
case where the values of the parameters are pushed in the direction that yields a large value 
of r]. We take A ^ 150 MeV, ttik ^ 4: MeV, 6m ~ -1.0 MeV, and C ^ 0.2. With 
these values we find that at T = 0.1 MeV (10^ K) the shear mean free path (Al) is 0.26 km, 
and r] = 1.7 x 10^"^ MeV^ = 2.3 x 10^^ergcm~^s~^. At this temperature the phonon's shear 
mean free path (see appendix A) is larger than the star, so the kaons provide the dominant 
contribution to the shear viscosity. 



hi{u) ^ 3.44 X lO-V^^ 
h2{iy) ^ 1.70 X 10-^z/ii 



(55) 
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FIG. 3: (Color online) The shear viscosity as a function of A. (See text for parameter values.) The 
points are calculated numerically. The straight lines are fits to the power law behaviors of (54). 
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FIG. 4: (Color online) The functions h\{y) and h2{v) (54). We scale out the u^^ power law behavior 
(see appendix C). On the rc-axis we show ly in units of the kaon velocity v ~ y^lSml/rriK in the 
non-kaon-condensed phase. 
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Finally, we check whether the regime of linear hydrodynamics is valid by evaluating the 
size of the corrections to the equilibrium stress-energy tensor. Linear hydrodynamics is 
appropriate if 5Tij <^ Tij (17). (Note that this is different from the criterion of validity for 
hydrodynamics in general, discussed in appendix A.) This inequality becomes 



where £ represents the length scale of the velocity gradients, and V is the typical fluid 
velocity which we assume is of order 1. If we use the energy density of free quark matter 
e 9/i^/(47r2) which is of order 10^° MeV^ at ^ 500 MeV, (and P < e, which is typically 
the case), and use the length scale £ ~ 1 km which is appropriate for oscillations of neutron 
stars, we find that linear hydrodynamics is valid as long as 



which is easily obeyed by the values of the shear viscosity that we have calculated for the 
CFL-KO phase. 

V. CONCLUSIONS 

In this paper, we have calculated the shear viscosity arising from self-interaction of the 
Goldstone kaon mode in the CFL-KO phase of quark matter. The shear viscosity from 
the other Goldstone mode, the superfluid phonon, has already been explored in Ref. [22]. 
We find the same temperature dependence that was found for the phonons in the 
CFL phase and for superfluid modes in a unitary Fermi gas [27]. Our final results are the 
approximate analytic expressions (54), (55) for the shear viscosity due to Goldstone kaons, 
and expressions (Al) and (A21) for the "shear mean free path" and "scattering mean free 
path" of the Goldstone kaons. 

Neutron star oscillations have a length scale in the kilometer range, so the phonon and 
Goldstone kaon fluids in a neutron star can only be described by hydrodynamics when their 
mean free paths are smaller than this. We argue that the shear mean free path is the 
appropriate quantity to use for this purpose (see appendix A). 

Because the coupling constants for the Goldstone kaons are roughly an order of magnitude 
larger than those for the superfluid phonon, the shear viscosity and mean free path of the 



ri < Vl{P + e) 



(56) 



7] < 10^^ MeV^ 



(57) 
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Goldstone kaon are both several orders of magnitude smaller than for the superfluid phonon 
(see Fig. 2). Using the shear mean free path (see appendix A), we find that the superfluid 
phonons in a neutron star are described by hydrodynamics at temperatures above about 
1 MeV (10^° K). The Goldstone kaons are hydrodynamic down to lower temperatures: the 
exact threshold depends sensitively on the value of the constants in the effective action, but 
could easily be lower than the 0.01 to 0.05 MeV range at which our treatment becomes 
invalid because the weak-interaction mass of the Goldstone kaon [21] must then be taken in 
to account. 

We conclude that, in the temperature range 0.01 MeV to 1 MeV, depending on the values 
of the coupling constants in their effective theory, Goldstone kaons may very well provide 
the dominant contribution to the shear viscosity in CFL-KO quark matter. 

There are several ways in which this work can be developed further. Firstly, we chose 
a specific form of the interaction lagrangian (15) which has one coupling constant, rather 
than the most general form (13) which has three; our calculation should be extended to 
the most general lagrangian. Secondly, it would be useful to extend our calculation to lower 
temperatures where, as noted above, one can no longer neglect the effects of weak interactions 
on the dispersion relation of the Goldstone kaon. It would also be interesting to study shear 
viscosity from light kaons in the non-kaon-condensed CFL phase. These particles were found 
to give a large contribution to the bulk viscosity even at temperatures as low as a tenth of 
their energy gap [1, 2]. Thirdly, we neglected scattering between the Goldstone kaons and 
the superfluid phonons. It would be interesting to see if these processes shorten the phonon 
mean free path and make a significant contribution to the shear viscosity. Fourthly, even 
though our calculation is open to extension and improvement in the ways just described, it 
would be interesting to perform an analysis along the lines of Ref. [34] to see whether the 
shear viscosity of the Goldstone kaons can have a significant effect on the development of 
r-modes in a quark star or hybrid neutron star. Fifthly, as discussed in appendix A, we did 
not consider how the interactions themselves would alter the dispersion relation. This could 
affect the calculation of the mean free path at leading order in the induced non-linearity, 
but would provide a subleading correction to the shear viscosity. Finally, even when the 
superfluid phonons or Goldstone kaons are not in the hydrodynamic regime, they can still 
transfer momentum over long distances, and it is important to investigate how they could 
provide ballistic-regime damping (as opposed to hydrodynamic viscous damping) of neutron 
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star oscillations. 
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Appendix A: Mean Free Path 

In this appendix we discuss the mean free path of the Goldstone kaon in the CFL-K° 
phase. We expect that hydrodynamics will be applicable to neutron star oscillations when 
the mean free path is well below the kilometer scale, since neutron star radii are about 
10 km. 

We study two definitions of the mean free path, which we call the "shear mean free 
path" l^^^^^ and the "scattering mean free path" Z^'^'**. The shear mean free path is based on 
the value of the shear viscosity itself, and is probably the physically relevant quantity for 
deciding when hydrodynamic calculations of shear viscosity are valid. The scattering mean 
free path is the average distance between collisions of the Goldstone kaons, including co- 
linear scattering events. Since co-linear scatterings do not contribute to the shear viscosity 
itself, it seems likely that this quantity is not the relevant one for finding the limits of validity 
of shear viscosity calculations. 



1. Shear mean free path 

We take our definition of the shear mean free path from Ref. [27], Eq. (48), where it is 
referred to as A^, 

^shcar _ _V 

n [p) 

where rj is the shear viscosity (given for kaons by (54) and (55)), (p) is the thermal average 
momentum, and n is the boson density, 

(p) = 2.7T/V , (A2) 

where u is the speed of the Goldstone bosons (given for kaons by (11)). 

This already allows us to make an estimate of the maximum shear viscosity that Goldstone 
kaons can provide, since it follows that t] ^ .3h'~'^T^l^^^^^ , so the maximum shear viscosity 
that could possibly occur in a neutron star at temperature T is when l^^'^^'^ ^ 1 km, i.e. 

'''""^ ^7 X 10-16 MeV ■ ^^^^ 
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For Goldstone kaons it is quite possible to get z/ ^ 0.1 by using small values of 6m. Using 
this value we find the following upper limits: at T = 0.01 MeV, ?7max ~ 10^^ MeV"^; at 
T = 0.1 MeV, r/„iax ~ 10^^ MeV^; at T = 1 MeV, r/^ax ~ 10^^ MeVl 

For superfluid phonons we can make more definite statements because there is less 
uncertainty about the parameters appearing in these expressions. The shear viscosity is 
?7 = 1.3 X 10~'^/ig/T^ and the speed u is generally assumed to take its perturbative value 
l/VS [22], in which case we immediately find that 



So phonon hydrodynamics becomes invalid in neutron stars when > 1 km, i.e. for 

/i ~ 500 MeV we require T > 1 MeV (10^° K), and at a temperature of 1 MeV the phonon 
shear viscosity is 5 x 10^'^ MeV^. 

2. Scattering mean free path 

Our calculation of this quantity closely follows that of Ref. [22]. The scattering mean 
free path is determined by the 2-body interaction cross section. The relevant scattering 
amplitude is the sum of four Feynman diagrams, the contact term and the s, t, and u- 
channel diagrams, see Fig. 1. In the s, t, and u-channels there is a virtual particle which 
can go on-shell, which means that its self-energy must be included to avoid an unphysical 
divergence. (The fact that the virtual particle can go on-shell means that the 2 — 2 collision 
rate already includes the contribution from 1 — )■ 2 splitting process, so these need not be 
calculated separately [35, 36]. We have performed this separate calculation and verified that 
the result has the same parametric dependence as the one we obtain below.) The rate (per 
unit volume) for the 2 — )■ 2 scattering process is 



'shear 
H 



4 X 10"^^ 



(A5) 




(2 nY6\P + K-P'- K') \M\' U fk (1 + fp') (1 + fk') 



(A6) 



where 




1 



(A7) 



and 



\M\' = \Ms\' + \Mt\'+\M^\' + I 



(A8) 
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where / represents the contact and interference terms, and the matrix elements are given in 
Eq. (48). The dominant contribution comes from |A^sP+|A^tP+|-^uP; since these each have 
a large enhancement when the virtual particle is close to being on-shell. The contact term 
and interference terms have no such enhancement, and so make a much smaller contribution 
to the rate. 

We now give a detailed explanation of the evaluation of |A^sP; the others can be obtained 
by similar methods. We first define the virtual particle momentum Q = {Eg, q) and shift 
the integral over k to an integral over q. We can then still use the momentum-conserving 
delta- function to do the integral over k'. We then choose the direction of q as the z-axis, 
so the angular part of the q integral gives a factor of 4 7r and the remaining integrand is 
azimuthally symmetric, so the 2 remaining azimuthal integrals give a factor of 27r. This 
leaves three integrals over the magnitudes of p , k' and q as well as over two polar angles, 
between q and p and between q and k'. We then introduce the auxiliary variable u via the 
identity 

/oo 
du5{u-Ep + Ep,)5{u + Ek-Ek'). (A9) 
'OO 

The integral over the two polar angles can then be done using these two delta- functions, 
leaving behind four integrals over p, k', u and q, 

POO POO p{u}+uq)/2u p(uj+vq)/2v 

Ts = A dq du dp dk' |A^,p Up f^^^p (1 + (1 + fuk') 

Jo J i> q J{u}—vq)/2u J{u]—uq)j2i> 

(AlO) 

To make the temperature dependence explicit we introduce a new set of variables (x, y, w, 
and z). For the s-channel these are 

up uk' u vq 

^ = Y y = ^ w = - ^^Y' ^ ' 

and the Mandelstam variable is s = uP' — z^. This leads to the full expression for the 
s-channel rate, 

^42^12 p .oo Pi^+.)/2 K^+z)/2 G{x,w,s)G{y,w,s) 

' dz / dw / dx dy r [x,y,w)- 



(A12) 

with 

F{x,y,w) = f^_^{l + fy){l + f^_y) G{x,w, s) = (3a;(l - u'^){w - x) - sf (A13) 
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The self-energy term has both a real and imaginary part, however the real part is much 
smaller [22]. We will therefore only consider the imaginary part Im n(cj,g), obtained in 
Refs. [1, 2, 22], 



Im n(w, z) = u+e{w^ - z^) + n-e(z2 - w;^^ 



,2 



1 



'w-\-z 
2 



g'^T^ 1 

n = ^ — TJI—r / dy G{y,w,t)fy[l + fy. 



where n"*" is relevant for the s-channel where w > z, and n~ is relevant for the t- and 
u-channels where w < z. We have neglected the tadpole contribution to the self-energy: it 
only corrects the kaon velocity by term proportional to AT^ [22]. 

In the remaining 4-dimensional integral (A12) we can now use a simple approximation to 
greatly simplify the integral. In the expression for Tg the integral over w is sharply peaked 
at the limit of integration where w = z, i.e. u = ±z/g (s = 0). The integral takes the form 

°° I{w,z) TlT^ 

^^(^2_^2)2 + (n+(u;,z)/T2)2 ^ 4m{z,z)z^^'^^ ' ^ ' 

This expression is valid when I{w, z) is slowly varying near the singular point w = z and 
when n^/T^ ^ 1 because 11 ~ g'^T^ and T ^ A,;^^'. Applying this approximation to 
the s-channel contribution, we find 

' 16z/i7(27r)5/.« 4 io 
where 



/ dz dx dyfj,-^{l + fy){l + f^^y) — ^ . (A15) 

Jo Jo Jo J--'- 



n+(^, = T^Klx^ I dx U fz-. z, 0) . (A16) 



We can then see that the integral over x in Eq. (A15) partially cancels the integral contained 
in n"*" leaving behind only a double integral. Using G{y,z,0) from Eq. (A13), we have 



- ' ' , ^ ^ > (AIT) 



where J is a pure number given by 



q r°° 

J = dz dyy'z^z-yfMl + fy){l + f,.y) 

^ 0.466 . (A18) 
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Applying the same method to the t- and u-channel integrals, we find that they all give the 
same rate, so the total rate is 



Ttotai = 1.40 .\ , ^ . (A19) 



The scattering mean free path is defined as 

7777 

where n is the particle density (A3). So the scattering mean free path of the Goldstone 
kaons is 

8 f4 8 ^2A4 

IT = 0.0881- = 3.44 X 10-^- — . , T'' , (A21) 

For comparison, the scattering mean free path of the phonon is [22] 

IT = 0.181 , = 5.02 X 10"=^^ (A22) 

so the ratio of the two is 

/scat 8 fi P\'^\()m\^ 

= 17.546 , i\ ^ 2.14 X 10-^ ^-, ' . (A23) 

Since /Xg is much larger than any of the other energy scales, this implies that the scattering 
mean free path of the kaon, like the shear mean free path, is generally much shorter than 
that of the phonon, giving the Goldstone kaon a much wider range of temperatures where 
it can be treated hydrodynamically. 

We have noted above that we expect the shear mean free path to be a better indicator 
of the range of validity of hydrodynamics, but for the sake of completeness we now estimate 
the temperature at which /|^^* will become greater than 1 km, in the case of very unfavorable 
parameter choices that lead to a long mean free path. We will use the values used at the 
end of Sec. IV B to illustrate how high the shear viscosity can be, namely = 150 MeV, 
A = 150 MeV, 5m = -1.0 MeV, tuk = 4.0 MeV and C = 0.2. In this case the scattering 
mean free path is shorter than 1 km for T > 0.006 MeV. For more favorable choices of 
the couplings this critical temperature will be much lower. In comparison, from (A22) the 
scattering mean free path for phonons is shorter than 1 km for T > 0.04 MeV. 
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We have used a linear dispersion for the kaon in calculating the mean free path, which 
is a requirement for getting a co-linear enhancement. However, there are sources of non- 
linearities in the dispersion. One comes directly from our expansion of the full kaon disper- 
sion in (11). If we had kept higher order terms, we would get a contribution that behaves 
as 



where 7 > 0. This positive curvature would still allow for the co-linear splitting and joining 
processes. Therefore, keeping this term would provide a subleading contribution to the 
calculation presented here. 

However, we ignored how the higher order derivative interactions themselves could change 
the kaon dispersion. Something similar has been calculated for the superfluid phonons, [37], 
where 7 was found to be negative and therefore the 1 o 2 processes are kinematically 
forbidden. If the corresponding non-linearity for the kaons were positive, then as above, 
the calculation presented here would remain the same. However, if the curvature were 
negative as for the phonons, then the mean free path would be altered at leading order. This 
is basically because the non-linearity itself would act to regulate the on-shell propagator 
and the scattering rate would go like I/7 instead of 1/H (where H is the self-energy). 
The appropriate scales to compare are 7T^ and H/T^ and in the case of the phonons, 
7T^ ^ n/T^, such that this correction would make the mean free path even larger and 
affect the validity of hydrodynamics. See [16] for a calculation involving the non-linear 
phonon dispersion and its affect on regulating the phonon propagator in a calculation of the 
thermal conductivity. However, it should be noted that including the non-linearity would 
only provide a subleading correction to the shear viscosity for either sign of 7 because the 
shear viscosity is insensitive to that region of phase space. 

Appendix B: Power counting sharply peaked integrals 

Here we discuss in more detail the evaluation of integrals of the type (A14), having a 
slowly-varying component / multiplied by a function with a sharp Lorentzian peak at the 
edge of the range of integration. In appendix A we assumed that / was non-zero at the 
edge of the range and we kept only the leading contribution. Here we include higher-order 




(A24) 
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corrections by Taylor-expanding the numerator, 




" I{z, z) + {w- z)I\z, z) + \{w- zfr\z, z) + ... 
(w"^ — z'^Y' + ^ 



(Bl) 



where l\z^ z) is the first derivative of / with respect to w, evaluated at w = z. We then find 



This gives us a scheme for power counting any integrals of the form given by (Bl). The 
relevant property is the dependence on e, since the collision integrals for transport properties 
take the form (Bl) with e = H/T^ oc g^T^ < 1. 

We can now justify the statement made in Sec. Ill that when we calculate the shear 
viscosity using a polynomial expansion of the function g{p) (38), the dominant contribution 
comes from g{p) = 1/p, i.e. choosing the minimum-exponent parameter n to be —1. 

Calculations of the mean free path and the shear viscosity both involve a rate calculation 
which contains collision integrals. In the mean free path collision integral (A12) there is 
a sharp peak in the integrand at w = z corresponding to a co-linear divergence, where 
two kaons have parallel momenta, and exchange a kaon whose momentum lies in the same 
direction. In the mean free path calculation this near-divergence is regulated by the self- 
energy, so the result depends on the self-energy (~ 1/e). 

In the case of the shear viscosity, we expect the integral not to have a co-linear divergence, 
since shear viscosity measures momentum transfer, so processes that do not change the 
momentum direction of the particles make no contribution. We therefore expect that the 
collision integrand in the shear viscosity should go to zero at w = 2; in such a way that 
the result does not depend on the self-energy. The true physical g{p) function will give 
an integrand that has this property. However, if we make a bad guess at g{p) (by using 
inappropriate basis polynomials in the expansion (38)) then each individual term will have 
a co-linear divergence, which will only cancel out when we add up the contributions from 
many terms. The best guesses for g{p) are therefore ones which yield a collision integrand 
with no co-linear divergence, i.e no dependence on the self-energy. 

In both the mean free path and the shear viscosity calculations the collision integral takes 
the form (Bl). The difference between them lies in behavior of the numerator I{w, z) in the 




(B2) 
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A/T i n - p X D o n p n t, 
parameter n 


Behavior of 
A?- near w = z 


^^-nPTiP'nnpnpp of 

collision integral 


n = —1 


{z — w)"^ 


independent 


n = -2 


{z — w) 


— In e 


n / -2,-1 


{z - wf 





TABLE III: Table of behavior of Ajj (part of the collision integrand) near the co-linear singularity, 
and the collision integral; e represents the self-energy. Only the n = —1 case has the proper physical 
suppression of co-linear contributions to the shear viscosity. 

co-linear regime w z. In the mean free path calculation (appendix A), the numerator 
stays finite in this regime, so the collision integral is of the form Ji (B2) and is strongly 
dependent on the self-energy. In the shear viscosity calculation, the numerator contains an 
additional factor A^^Ajj (see Eq. (43)), which we'll call the transport term. The behavior of 
the transport term in the co-linear regime is therefore crucial in suppressing the self-energy 
dependence. 

Using expansions of g{p) with only one term (A^ = 0) we calculate A?- and the shear- 
viscosity collision integral for different choices of n. We summarize the results in Table III. 
We see that the choice n = —1 fully suppresses the co-linear singularity and gives a collision 
integral that is independent of the self-energy. The choice n = —2 partially suppresses 
co-linear scattering and gives the collision integral a very weak dependence on the self- 
energy. Other values of n do not suppress the co-linear scattering at all and are akin to 
the calculation of the mean free path. This explains our finding in Sec. Ill, Table II that 
n = — 1 is the the optimal choice for fast convergence of the polynomial approximation for 
g{p) in the shear viscosity, that n = —2 is the next best choice, and other values of n have 
very poor convergence. 

Appendix C: Approximate evaluation of the collision integral 

Here we describe how the collision integral is reduced to a five- dimensional numerical 
integral, in which we have factored out the temperature dependence and part of the de- 
pendence on the kaon speed u. We begin with the matrix M (43) that enters in to the 
calculation of the shear viscosity. As described in Sec. Ill and appendices A and B, we get a 
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good estimate of the shear viscosity by assuming g{p) = 1/p, i.e. we set = and n = —1. 
Then, as described after Eq. (49), we can ehminate seven of these integrals by using the 
(5-function and spherical symmetry. We rescale the momenta with temperature, and find 

Moo = (^) / /. /. (1 + /.) (1 + U) IM^, 9, A) P A°. A°. (CI) 

where fx = l/{e^ — l), and 



Jo J~i J -I 



1 /■! ^2 



dV = dx dy da dl3 d(j) , (C2) 



a 



pp fk vk' xyil—a) up' 

X = — V = — z = — — = — — ; w = —— = X + y — z , (C3) 

T y 7- T x(l - /3) + ?/(l - 7) T ^ ' V ; 



a = p.k /3 = p-k' 7 = k-k' = a/3 + V^(r^l^2y(Y^r^cos(0) (C4) 

From (48), Ai = {h'/T)^Ai depends on the speed u and the couplings g and A as well as 
the rescaled momenta x, y, z. From (44), = (z//T) A^^ depends only on the rescaled 
momenta. The expression for z comes from solving the energy-conserving (5-function, and (p 
is the difference in azimuthal angles between the vectors p, k and p, k'. 

We have scaled out all the temperature dependence of the integrand, but there is still 
some dependence on u and the couplings g and A which comes in via Ai. The integral can 
be evaluated numerically using (44) and (48) for given values of g, A. 

We now show how to obtain the approximate analytic forms for 77 given in (54), which 
are valid in the regime g'^/X <^ 1 and g'^/X ^ 1. 

Using (CI) in (46) and evaluating Aq^ from (39), we find that 

const 

where const represents a function independent of u. We can obtain the function /ii(i^) in (54) 
by going to large A in which case Ai = M.ci and then calculating the shear viscosity with 
all the dimensionful parameters in the coupling constant A set equal to unity. To calculate 
h2{u), we go to large g where JH = Ais + Ait + -Mu and do the same thing. In both cases 
we find 

^^'^ - JdTm\M,,iuW ^^'^ 
where Ai^2 is a pure number and f(T) represents the parts of the integrand of (CI) that are 

independent of u. From (43) we find that 

-Mi,2 ~ cfi + cS>' + cfy (C7) 
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Therefore, we expect 



,11 



y-4 (i) 2i ' 
l^i=0 "1,2*^ 



(C8) 



justifying our statement in the paragraph below Eq. (54) and the resulting scaling in Fig. 4. 

Finally, we can explain why hi and /i2 are so small (see Fig.4). This is a direct result of Mqo 
being large. Because the all interactions of the Goldstone kaons are derivative interactions, 
the collision integral involves high powers of momenta. Schematically, it has the form 



Jo 

where c? = 12 in our case. Since 1/13! ~ 1.6 x 10^^" it is not surprising that hi and h2 are 
of that order. 
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